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Ph. ABSTRACT 

The high redshift radiogalaxy 1243+036 (z = 3.6) presents an asymmetric Lya 
fSj ; profile of FWHM 1550 kms -1 as measured by van Ojik et al. We propose that the 

blue asymmetry in the Lya profile is not due to narrow absorption dips but consists 
of narrow emission peaks. We interpret the blueshifted peaks near —1130, —850 and 
—550 km s" 1 (relative to the peak of full profile) as being the result of Fermi acceleration 
of Lya produced by jet-induced star formation in the wake of a 300 kms -1 shock. This 
■^j- ■ shock would be caused by the deflection of the radio-jet at the observed position of 

the radio bend which also coincides spatially with the excess Lya emission reported by 
\ van Ojik et al. 
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1. Introduction 



Lyman a is the strongest line observed in 
very high redshift radio galaxies (HZRG). Al- 
though the brightness of Lya peaks near or at 
the nuclear position, most of the emission is 
spatially resolved with the fainter emission ex- 
tending up to radii 40-130 kpc. The Lya pro- 
file in HZRG is characterized by a FWHM in 
the range 700-1600 km s _1 (van Ojik 1995 and 
references therein). The intermediate resolution 
study of Lya profiles carried out by van Ojik 
et al. (1997) has revealed the presence of troughs 
which are well explained by Hi gas absorption 
present in the environment of the parent radio 
galaxy. In their sample, however, the HZRG 
1243+036 (z = 3.6) may call for a different 
interpretation, namely that the broad profile 
presents true narrow emission features in be- 
tween the narrow 'dips' which van Ojik et al. 
(1996: v096) have interpreted as absorption fea- 
tures. The repeated scattering of the resonance 
Lya line across a shock discontinuity was shown 
by Neufeld & McKee (1988: NM88) to result in a 
systematic blueshift of Lya. In this paper, we de- 
velop in more detail the Fermi acceleration model 
and propose that the narrow features observed 
by v096 correspond to a small number of across- 
shock scatterings. 

2. The model 

The calculations of the Lya profile is carried 
out in two steps. First we compute the emergent 
profile from the ionized gas, second we consider 
the effect of Fermi acceleration across a shock 
discontinuity which is taking place between the 
observer and the Lya emission gas. 

2.1. Double-horned source profile 

A resonant photon escaping from a static slab 
will diffuse in frequency and space until it can 
escape in the wings of the opacity profile. This 
leads to a symmetric double-horned profile with 
a separation, v sep , between the two peaks of: 



v sep « 390 (A^V^) km s_1 . (!) 

where 10 21 A^|i u cm~ 2 is the neutral hydrogen 
column density surrounding the emitting ionized 
gas and 10 2 T.f ou K its temperature. We properly 
take into account the skewness at large iV|° u of 
the two emission peaks (cf. Bonilha et al. 1979) 
by using the numerical prescription described in 
Appendix B of Binette et al. (1993). We can ei- 
ther assume the ionized gas to be homogeneous 
and isothermal (structure-less) or we can inte- 
grate the Lya profile escaping a cooling shock 
slab or a stratified photoionized slab. In either 
case, we assume a plane-parallel geometry for the 
transfer which allows two possible perspectives 
for an external observer: front or back. The 
front perspective corresponds to the precursor 
side (shock wave case) or the UV source side 
(photoionization case). The fraction of Lya es- 
caping towards the front or the back can be quite 
different, especially for an ionization bounded 
photoionized slab or for a shock wave with a trail- 
ing thick shell. The H I gas beyond the Stromgren 
boundary or within the trailing shell acts like 
mirrors which favor Lya photons to escape over- 
whelmingly towards the front. This effect as well 
as Lya destruction due to internal dust are fully 
taken into account using the formalism developed 
in Appendix A of Binette et al. (1993). We label 
iV™ r the column density of the trailing H I shell. 

2.2. Fermi acceleration of resonant pho- 
tons 

As for the effect of Fermi acceleration on the 
Lya profile, unless stated otherwise we will con- 
sider the front perspective whereby the shock dis- 
continuity is taking place between the observer 
and the region producing the line. This is dif- 
ferent from the (back) perspective^ in NM88 in 

Our front perspective prevents the components (a) and 
(b) discussed in NM88 from being visible. 
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which the shock was on the far-side of the emer- 
gent Lya (see their Fig. 1). As in NM88, we treat 
the slabs of scattering material on either side of 
the shock front Nl^ and iV™ T ' as two converg- 
ing, partially reflecting mirrors. Each passage of 
a resonant photon across the shock front corre- 
sponds to a reflection and, for an isotropic radi- 
ation field, the mean blueshift after n passages 
across the discontinuity is nV s /2, where V, is the 
shock wave velocity. If the transmissivity of the 
mirrors on the kth. reflection is T k and the reflec- 
tivity is 1 — Tfc , the fraction of resonant photons 
escaping at exactly n shock crossings is 

f{n) = T n l[ (1 - T k ) « T n exp(- f T k dk) . 

(2) 

The right term is a valid approximation only in 
the limit of large n and < 1. Unlike NM88, we 
have used instead the middle term in which each 
crossing of the discontinuity is calculated follow- 
ing the algorithm described in Appendix A. The 
main advantage is that it allows us to cover the 
case of a small number of reflections. Further- 
more, it retains higher resolution of the spectral 
features appearing on the transmitted profile and 
it allows a more accurate treatment of the effect 
of dust absorption (see Appendix A). We define 
fj,Y) as the dust-to-gas ratio relative to the local 
interstellar value. 

The profile of non-resonant lines of species i 
of atomic mass mi are assumed gaussian with a 
FWHM given by 1.665 (2kT/rrii + v% urb ) 1/2 where 
vturb is the local 'most probable' turbulent ve- 
locity. This summarizes the implementation of 
profile determination in the multipurpose code 
mappings IC (cf. Ferruit et al. 1997). 

2.3. Profiles produced by Fermi acceler- 
ation 

To illustrate the type of profiles obtained, we 
have made calculations under a wide range of in- 
put parameters. In these calculations, the shock 



velocity, 100 Vf kms , is constant at V-f = 3.0 . 
Increasing Vf does not necessarily change the 
shape of profile if the opacities are increased ac- 
cordingly (although it does blueshift the profile 
further). For definiteness, we have set = 
N™ 1 " in all our calculations^]. We adopted a tem- 
perature of 100 K for the gas responsible for the 
scattering. The choice of temperature impacts 
only on v sep [see equation (1)] and therefore only 
the product N^y/lf^ matters to the calcula- 
tions. At constant temperature, the two essential 
free parameters are therefore A^™^ (= N^ ) and 
A^2i" [the column depth immediately surround- 
ing the ionized g equation (1)]. The dif- 
ferent profiles resulting from Fermi acceleration 
are presented in Fig. 1 and the corresponding in- 
put parameters in Table 1. It is seen from Fig. 1 
that the higher N™ 7 ', the higher the blueshift 
of the emergent Fermi profile (see Appendix B 
for a comparison with NM88). v sep is an im- 
portant factor governing the appearance of the 
profile when it is ;> V s . The effect of dust on 
Fermi acceleration is not as severe as in NM88. 
The main reason is that we considered only the 
effects of dust within the two mirrors and its im- 
plication on the Fermi profile per se while dis- 
regarding the presence of dust within iVIf", the 
gas immediately surrounding the source. In any 
event, little dust absorption would take place 
when N2° u is small, a possibility strongly favored 
in the discussion below. Also we expect the ex- 
tremely hot postshock gas surrounding the em- 
bedded condensations to be devoid of dust. In 
our calculations of the dust effects on the Fermi 
profile, we have taken into account the mean ge- 
ometrical depth traversed at each shock crossing 
by the resonant photon. Fig. lj illustrates the 
case with /io = 0.2 (dust content of mirrors only) 
which should be compared with the dustfree case 
of Fig. li. 



6 In the event of a large imbalance, clearly Lya would escape 
in the direction of the thinnest of the two mirrors. 
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3. Analysis of the 1243+036 data 



4. Fermi acceleration in 1243+036 ? 



We have reduced the ESO archived observa- 
tions of PKS B1243+036 taken at the NTT with 
EMMI by v096 on 14 April 1994. The 4 one 
hour duration long-slit spectra with PA=152° 
(aligned with the radio axis) were registered on 
a 2048 x 2048 Tektronix CCD. The 2x2 binning 
gave a spatial scale of 0.54"pix _1 . Grating 6 was 
used with a wide slit of 2.5". Our analysis of 
the arc spectra indicates a resolution (FWHM) 
of ~ 2.6A. 

After performing a flat-field division, sky sub- 
traction, flux calibration and correction for at- 
mospheric extinction using the IRAF reduction 
package, we extracted a one-dimensional spec- 
trum using a restricted window of 7", centered on 
the bending point of the radio source (cf. Fig. 4 
and 14 of v096). The resulting spectrum is plot- 
ted as solid line in the upper frame of Fig. 2. It 
differs from that of v096 result of consider- 
ing a window which is both smaller and centered 
around the radio bend. 

Several small peaks (present in the individual 
exposures) appear on the Lya profile at 5536, 
5541 and 5547A. In order to isolate these fea- 
tures, we smoothed the spectrum and subtracted 
the result from the original spectrum. The re- 
sult is plotted as a dotted line in the upper 
frame of Fig. 2. Obviously, the negative part of 
this curve is unphysical whether the features are 
emission peaks or absorption dips. To remedy 
this, we supposed that the underlying Lya pro- 
file is intrinsically symmetric since almost all the 
Lya profiles measured by van Ojik et al. (1997) 
turned out to be symmetric after correction of 
the absorption troughs. Taking the red side of 
the smoothed profile as our reference, we folded it 
around the peak position (\m<ix = 5557.2A) over 
the blue side and subtracted it from the original 
profile. The resulting spectrum is shown in the 
lower frame of Fig. 2 (dotted line) . 



Are the high frequency features (dotted line 
in Fig. 2b) absorption dips or emission peaks? 
The absorption interpretation is ambiguous be- 
cause the V-shaped dips are unresolved in wave- 
length, i.e., their widths are less than the width 
of the flat-topped instrumental profile. This can 
be seen in Fig. 2b where the instrumental profile 
is shown in the inset box. On the other hand, 
if interpreted as emission features, these features 
appear to fit well the instrumental profile, as is 
evident from Fig. 2b. Furthermore, the peaks 
are equally spaced in wavelength and consider- 
ably blue shifted relative to systemic velocity, all 
in accordance with Fermi acceleration. Clearly, 
higher spatial and spectral resolution spectra 
would help settle the question. The existence of 
a shock is not in doubt, however, since 1243+036 
clearly shows the presence of a radio bend at 
a distance of 2" from the nucleus (13.6 kpc for 
Ho = 50 km s _1 Mpc _1 and qo = 0.5), that is 
at the same spatial position where v096 reports 
genuine enhanced Lya emission (underlying their 
absorption dips). 

Let us thereafter assume that the narrow emis- 
sion peaks are independent from the underlying 
symmetric Lya profile and that they are emit- 
ted near or at the radio bend. Between suc- 
cessive peaks, two reflections must necessarily 
take place. Therefore the velocity of the shock 
in our model is fixed by the distance between 
peaks which gives Vf ~ 3.0 . If the emission 
peaks corresponded to gas photoionized by the 
nucleus, we could not explain these by Fermi ac- 
celeration since one of the necessary H I mirrors 
(at the back) could not exist between the ion- 
ized gas and the source. Our Fermi acceleration 
model therefore requires in situ photoionization, 
either by the shock itself (a) or by an indepen- 
dent local source (b). 

In case a, the post-shock temperature is of 
order 1.5 x 10 6 K and in order to reproduce 
the integrated luminosity of the emission features 
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(2.4 x 10 43 erg s^ 1 ), one requires the uncomfort- 
ably high pre-shock density of no = 3 x 10 6 cm -3 
for a shock front large as 10 4 pc 2 (we in- 

cluded the ionizing photons generated within the 
shocked gas). Furthermore, the abundances of 
the heavy elements must be extremely low since 
C IV is absent from the spectra (v096) . When we 
compute the Lya profile from a shock, it turns 
out significantly wider (after convolution) than 
the observed narrow peaks. 

In case b, Lya emission originates from star- 
bursts triggered by the passing of the shock wave 
which is the result of the deflection of the radio 
jet. This possibility goes along the suggestion 
of Rees (1989) that enhanced star formation can 
be triggered by the interaction of the radio jet 
with the environment of HZRG. By inspection of 
many profiles such as in Fig. 1, two possibilities 
arise. One is that the Hn gas is embedded in 
thick neutral gas such that v sep is equal to V s . In 
this particular case, the convolved emission peaks 
are sufficiently narrow. Such a model is shown in 
Fig. 11 and it gives, after convolution, an excel- 
lent fit. A less ad hoc situation which does not 
rely on having v sep = V s is to have the H n regions 
not fully embedded but rather in their cham- 
pagne phase (Tenorio-Tagle 1979). This geom- 
etry which facilitates considerably the escape of 
Lya into the surrounding shocked gas results in a 
small v sep . Such a nebular model was computed 
using a stellar atmosphere model of 40 000 K. The 
Fermi profile obtained gives an encouraging fit 
as shown by the solid line in Fig. 2b. It was 
computed using Vf = 3.0 and N%± r = 1.3 and 
has a very similar profile to the simpler model of 
Fig. Ik when the latter is convolved. A point in 
favor of the Fermi interpretation is that the peaks 
appear at —4, —3, —2 V s relative to the adopted 
value of line center at 1.3A (70kms _1 ~ jV s ) to 
the red of Xmax , while absorption dips ought to 
be distributed randomly in velocity space. 

It is interesting to note that the shock crossing 
time across a projected distance of 0.2" -an up- 
per limit to the projected radio hot spot size- is 



~ 4 x 10 6 years. Incidentally, a preshock gas den- 
sity of 0.16 cm -3 [eight times higher than in Rees 
1989] with a column of N%[ w = 1.3 possess a geo- 
metrical depth of 0.2". After being shocked, such 
gas would cool on a timescale of ~ 0.8 x 10 6 years 
from a postshock temperature of 1.5 x 10 6 K. If 
the flow became very chaotic after the gas has 
cooled, this might preserve a small value for v se p 
even if the H n regions were to turn on at a later 
stage further downstream. 

Interestingly, the only other candidate HZRG 
for which Fermi acceleration was proposed, 3C326.1 
(NM88), is similarly devoid of Civ emission, a 
fact difficult to reconcile with either nuclear pho- 
toionization by a powerlaw or with high velocity 
(photoionizing) shocks. Jet-triggered star forma- 
tion as we argued for above would solve this co- 
nundrum since C iv is then expected to be very 
weak even for solar metallicity. 

We wish to thank Huub Rottgering for direct- 
ing us to the archived NTT data of 1243+036 
which were reduced a second time by one of 
us (BJ). We also thank Irapuan Rodrigues de 
Oliveira for an initial quick survey on the type 
of profiles produced by Fermi acceleration. One 
of us (JCLW) acknowledges support by NASA 
Astrophysics Theory Program grant NAG5-3836 
and by a Graduate Research Board award from 
the University of Maryland. 

A. Numerical algorithm to calculate Fermi 
profiles 

We discretize all profiles in vectors whose in- 
dices are proportional to the velocity shift (vi = 
i 5v) from line center Ao- 6v is a proportional- 
ity constant sufficiently small to allow fair sam- 
pling of the source profile. The source line profile 
emergent at the shock discontinuity is defined as 
<p(vi) (cf. § p.l|) . It is centered and symmetric 
around r>i = 0. In this scheme, the observer is ei- 
ther on the side of the precursor ('front' perspec- 
tive): p = +1, or downstream beyond the trail- 
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ing Hi shell ('back' perspective): p = — 1. After 
j + 1 passage across the discontinuity of the pho- 
tons coming from position Vi in the source profile, 
their contribution to the observed profile ^(vk) 
for an observer whose perspective is p = (— l)i is 
given by: 

¥>(OT(^)n (i-r(«n)) , (Al) 

n=0 

where A: is the nearest integer to Vj/8v, and with 
-Uj and v n given by: 

fj = «i-jVi/2-(-l) J V 8 , (A2) 

v n = Vl -nV s /2-(l + {-iy)V s /2 . (A3) 

The profiles (p +1 and c^ -1 are obtained by adding 
contributions from all j until 1 — T(vj) becomes 
negligible (note that j starts at zero). The 
transmissivities are approximated as T(v n ) = 
[1 + 0.75r(v ra )] _1 . Assuming a Holtsmark func- 
tion (Lang 1974), T(v n ) is the line opacity away 
from line center (by the amount \ov n /c) due to 
the opacity of either the precursor (n even) or 
of the shell (n odd). The line center opacities 
To are derived from N^ (if p = +1) and iV™ r 
(if p = —1). The attenuation of the transmissiv- 
ity due to dust absorption has been taken into 
account by determining the column density of 
dust crossed by the resonant photon at each pas- 
sage across the discontinuity (assuming that re- 
flection takes place at a geometrical depth such 
that r(v n ) « 1). As in NM88, we neglect the dif- 
fusion in frequency that results from scattering 
within the static material on either side of the 
shock front. 

B. Comparison with NM88 

The position of the peak intensity within the 
Fermi accelerated profile, V peak , shifts towards 
the blue as N^ 1 " and/or Vf are progressively 
increased. Our numerical calculations lead to 
the following relation valid only when N™ 1 " » 
0.002(U/) 2 : 



ypeak _ ^ ^ + ^ (ys N ^ I ? ^ 

with 5 defined as: 

5 = 0.16 (iV£n-3 y/vf , (B2) 

where 100 vf eak kms -1 is the peak intensity ve- 
locity. This expression assumes that = 
N%l ir . When N^ ir £ 10 4 (that is, for a negli- 
gible 5), the above expression for V 7 peah becomes 
identical to that of NM88. A more general ex- 
pression is to appear in Joguet & Binette (1998). 
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Table 1 

Parameters characterizing profiles <p +l , ip^ 1 of Fig. 1. 



Model 


Vf 


ATSOU 

IV 21 


ATmir a 
JV 21 


n^sou 


rpmir 
1 2 


MD b 


a 


3.0 


0.01 


0.1 


1.0 


1.0 





b 


3.0 


0.2 


0.2 


1.0 


1.0 





c 


3.0 


0.1 


1.0 


1.0 


1.0 





d 


3.0 


1.0 


1.0 


1.0 


1.0 





e 


3.0 


5.0 


1.0 


1.0 


1.0 





f 


3.0 


1.0 


0.1 


1.0 


1.0 





g 


3.0 


0.1 


0.01 


1.0 


1.0 





h 


3.0 


0.1 


0.001 


1.0 


1.0 





i 


3.0 


3.0 


1.0 


1.0 


1.0 





j 


3.0 


3.0 


1.0 


1.0 


1.0 


0.2 C 


k 


3.0 


0.0002 


1.0 


100 


1.0 





1 


3.0 


0.28 


0.53 


1.0 


1.0 






a We have set Nl^ ec = NJ^ r in all the above calculations. 

b The dust-to-gas ratio within the two mirrors relative to the 
local interstellar value (p,-o=l). 

b The fraction of Lyoj absorbed by dust within the two mir- 
rors is 63%. 
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Fig. 1. — The dash-dotted line represents dif- 
ferent double-horned source profiles (cf. § |2,1[ ) 
as 9j function of velocity from line center (v = 
c (A — Ao)/Ao). The profiles resulting from Fermi 
acceleration were calculated using the parame- 
ters given in Table 1. After multiplication by 
an appropriate constant (e.g. X4, to provide bet- 
ter comparison with the source profile), they are 
plotted according to the observer's perspective: 
front (tp +1 )~- solid line or back (ip~ 1 ): dotted line. 
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Fig. 2. — Top: The observed Lya profile of 
1243+036 (solid line). The dotted line repre- 
sents the high frequency structures extracted by 
smoothing and subtraction. Bottom: The differ- 
ence between the left half and the smoothed and 
folded right half of the Lya profile is shown as 
dotted line. The lower ordinate is the velocity 
shift if we adopt as line center Ao = 5558. 5A. 
The solid line represents the Fermi acceleration 
of the Lya profile of an Hn region in its cham- 
pagne phase. This profile has been convolved by 
the instrumental profile shown in the inset boxes. 
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